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Abstract 

The Monte Carlo (MC) Hamiltonian is a new stochastic method to solve many- 
body problems. The MC Hamiltonian represents an effective Hamiltonian in a fi- 
nite energy window. We construct it from the classical action via Monte Carlo 
with importance sampling. The MC Hamiltonian yields the energy spectrum and 
corresponding wave functions in a low energy window. This allows to compute ther- 
modynamical observables in a low temperature window. We show the working of 
the MC Hamiltonian by an example from lattice field theory (Klein-Gordon model). 

Key words: Many-body systems, thermodynamical observables, Monte Carlo 
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1 Introduction 

Many interesting phenomena in physics occur in many-body systems. The 
theoretical difficulty lies in solving models of many-body systems. A brilliant 
idea is the renormalisation group approach a la Kadanoff and Wilson [1]. It 
suggests to compute critical phenomena by constructing an effective (renor- 
malized) Hamiltonian, by thinning out degrees of freedom in an appropriate 
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manner. Here we suggest to adapt this idea for the purpose to construct an ef- 
fective Hamiltonian in a low-energy window. The basic idea to thin out degrees 
of freedom is taken over from lattice field theory, where infinite-dimensional 
integrals (path integrals) are successfully simulated by a small number of rep- 
resentative configurations. In analogy to the representative configurations of 
path integrals, we construct a representative basis of Hilbert states via a Monte 
Carlo (MC) algorithm with importance sampling. Furthermore we use the MC 
algorithm to compute matrix elements of the transition amplitude. 

In recent years, Monte Carlo methods have been widely used to solve problems 
in quantum physics. For example, with quantum Monte Carlo there has been 
improvement in nuclear shell model calculations [2]. A proposal to solve the 
sign problem in Monte Carlo Greens function method, useful for spin models 
has been made by Sorello [3]. Lee et al. [4] have suggested a method to di- 
agonalize Hamiltonians, via a random search of basis vectors having a large 
overlap with low-energy eigenstates. In contrast to that, in this work we con- 
struct matrix elements of the transition amplitude and hence the Hamiltonian 
via the path integral starting from the action. 



2 Monte Carlo Hamiltonian 

Let us discuss the construction of the MC Hamiltonian introduced in Ref. [5] 
in several steps. First, consider in 1-dim quantum mechanics the motion of 
a single particle of mass m under the influence of a local potential V(x). Its 
classical action is given by 



Given the classical action, one can determine quantum mechanical (Q.M.) 
transition amplitudes. Similar to the approach in Lagrangian lattice field the- 
ory, we use imaginary time in what follows. We consider the transition am- 
plitude for some finite time T (t in — 0, tfi — T) and for all combinations of 
positions Xi, Xj. Here Xi, Xj run over a finite discrete set of points {xi, . . . , xn} 
located on the real axis. Suppose these points are equidistantly distributed 
(spacing Ax) over some interval [-L, +L\. The transition amplitudes are given 
by the (Euclidean) path integral, 



where Se denotes the Euclidean action. Suppose for the moment that those 
transition amplitudes were known. From that one can construct an approxi- 
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mate, i.e. effective Hamiltonian. Note that the matrix M(T) = [Mjj(T)]jv x jv 
is a real, positive and Hermitian matrix. It can be factorized into a unitary 
matrix U and a real diagonal matrix D(T), 

M(T) = C/t D(T) U . (3) 



The matrix M(T) can be expressed in terms of the full Hamiltonian H by 

M l3 {T) = { Xi \e- HT / h \ Xj ) . (4) 

The matrices U and D have the following meaning, 



Ul=< Xl \E e k ff > , 

D k (T) = e- E l SfT l h . (5) 

The eigenvector > of the effective Hamiltonian H e ff can be identified 

with the column k of matrix W and the energy eigenvalues E^f* of H e ff can 
be identified with the logarithm of the diagonal matrix elements of D(T). 
Thus via the factorisation (3) one can define an effective Hamiltonian, 

N 

H eff = Y,\E e k ff >E e k ff <E^\. (6) 
k=i 



Note that in the above we have been mathematically a bit sloppy. The states 
\xi) are not Hilbert states. We have to replace by some "localized" Hilbert 
state. This can be done by introducing box states. We associate to each X{ some 
box state bi, defined by 



h(x) = < 



l/\fKxl if Xi < x < x i+ i ^ 
else 



where Axi = — X{. When we use an equidistant distribution of Xj, i.e. 
Axi = Ax, we refer to the basis of box states as the regular basis. Below we 
will consider a stochastic basis. 



3 Matrix elements 



We compute the matrix element Mjj(T) directly from the action via Monte 
Carlo with importance sampling. This is done by writing each transition ma- 
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trix element as a ratio of two path integrals. This is done by splitting the 
action into a non- interacting part and an interacting part, 

S = So + S v . (8) 



This allows to express the transition amplitude (2) by 

(0) J[dx] exp[-S v [x]/h] exp[-S [x]/h}\l" T 
Mij(T) = M§\T) ^— ! l T J '° . (9) 



f[dx] exp[-S Q [x]/h]\ 



Xj,0 



Here e~ s °/ h is the weight factor and e~ Sv l h is the observable. (T) stands 
for the transition amplitude of the noninteracting system, which is known ana- 
lytically. For details see ref.[5]. Carrying out these steps allows us to construct 
an effective Hamiltonian, which reproduces well low energy physics, provided 
that the nodes x^ cover a large enough interval (depending on the range of the 
potential) and the resolution Ax is small enough. This can be achieved with 
a small numerical effort for a one-body system in 1 dimension. Our goal is, 
however, to solve many-body systems. What to do in such case is the subject 
of the following section. 



4 Stochastic basis 



It is evident that the regular basis defined above becomes prohibitively large 
when applied to a many-body system. For example, in a spin model of a 
1-dimensional chain of 30 atoms with spin 1/2, the Hilbert space has the 
dimension D = 2 30 = 1073741824. For such situations we wish to construct a 
smaller basis which gives an effective Hamiltonian reproducing well low-energy 
observables. Why should such a basis exist in the first place? The heuristic 
argument is that the Euclidean path integral, when evaluated via Monte Carlo 
with importance sampling, gives a good answer for the transition amplitude. 
In particular, this is possible by taking into account a "small" number of 
configurations (e.g. on the order of N stoch = 100 — 1000). In a crude way the 
configurations correspond to basis functions. Thus we expect that suitably 
chosen basis functions exist, the number of which is in the order of 100 - 1000, 
which yields a satisfactory effective low energy Hamiltonian. 

We shall construct a small basis in a stochastic way, using Monte Carlo. We 
proceed by the following steps: (i) Compute the Euclidean Green's function 
Ge(x, T; 0, 0). We define P(x) as the probability density (normalized to unity) 
obtained from Ge(x, T; 0, 0). (ii) Find an algorithm giving a random variable x 
distributed according to P(x) and draw samples from this distribution, giving 
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nodes, say x v . Finally, one obtains the stochastic basis by constructing the 
corresponding box states from the nodes x v . This goal can be achieved in an 
elegant and efficient manner via the Euclidean path integral, expressing P(x) 
as 



S[dy]e^[-S E [y]/hf' 
P(x) = 



0,0 



S3* ■ (io) 

J^dxJ[dy}exp[-S E [y}/h] 

0,0 

Using a Monte Carlo algorithm with importance sampling (e.g., Metropolis [6]) 
one generates representative paths, which all start at x = 0, t = and arrive at 
some position x at time t = T. Let us denote those paths (configurations) by 
C v = x v {t). We denote the endpoint of path C u at time t = T by x v = x v (T). 
Those form the stochastically selected nodes, which define the stochastic basis. 
One should note that the Green's function Ge, the probability density P and 
hence the stochastic basis all depend on the choice of the parameter T, related 
to the inverse temperature (3 via (3 = T/h. It turns out that T influences 
the size of the finite temperature window in thermodynamical observables. 
However, it is not the only parameter to do so, as this window also depends 
on the size of the basis. 



5 Numerical results 



5.1 Quantum mechanics 



The Monte Carlo Hamiltonian has been tested on a number of quantum me- 
chanical potential models, by computing the spectrum, wave functions and 
thermodynamical observables in some low temperature window [5,7-9]. Al- 
though the results from the Monte Carlo Hamiltonian agree well with the ex- 
act results, low-dimensional models in Q.M. are not the target of this method. 
The target are rather high-dimensional models or models with a large number 
of degrees of freedom. Our expectation that the MC Hamiltonian works well 
in such system is guided by the analogy of numerical integration with Monte 
Carlo: The Monte Carlo method is not the most efficient method when doing 
low-dimensional integrals. However, it is most efficient when doing integrals 
in high dimensions (e.g. path integrals). 
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5.2 Klein-Gordon model 



We consider in D=l a chain of coupled harmonic oscillators, which is equiva- 
lent to the Klein-Gordon field on a 1 + 1 lattice. The model is given by 



S = J dtT-V 

n=l Z 

V= l -Y.n 2 {q n -q n+l ) 2 + nlq 2 n . (11) 
Z n=l 

The parameters have been chosen as m = 1, f] = 1, Sl = 2, N osc = 9, a = 1, 
T = 2 and h = 1. After the stochastic basis with N stoc h is generated, we 
obtain the matrix elements M n / n Then we compute the eigenvalues and eigen- 
vectors using the method described in Sect. (2). Table [1] gives a comparison 
between the spectrum from the effective Hamiltonian with the stochastic ba- 
sis and the analytic result for the first 20 states. We display the results for 
Ngtoch = 10, 100, 1000 and observe good agreement for large enough N stoc h- 
This means that the stochastic basis works well. 



We have also computed thermodynamical quantities such as the partition 
function Z, free energy F, average energy U = E and specific heat C. The 
analytical results are 



Z(f3) = Tr(exp(-(3H))= JJ 



Nose 



1-- 



\ 2sinh (Phcui/2) ' 



E{p) = -Tr (H exp {-(3H)) = = £ ^ coth (/3^/2) , 



= or = ~ k BP m=^L ^XT^T 7^ > ( 12 ) 



^ c / ^/2 x 2 

where 



^ = ^fig + 4i2 2 sin 2 (p/Ax/2) , 
Ap = 2tt/(A^ osc Ax) , 
x j = [-(N OBC -l)/2 + (j-l)]Ax , 

Pi = [-(N 0BC -l)/2 + (Z-l)]Ap. (13) 
Here j and / run from 1 to iV osc (number of oscillators). Ax = a = 1 is the 
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Table 1 

Spectrum of Klein-Gordon model on the lattice, MC Hamiltonian (stochastic basis) 
versus exact result. 



n 


Ef (Nstoch = 10) 


Ef (N stoch = 100) 


Ef (N stoch = 1000) 


^exact 


1 


10.671 


10.960 


10.905 


10.944 


2 


12.725 


12.964 


12.957 


12.944 


3 


13.105 


13.109 


12.985 


13.057 


4 


13.186 


13.198 


13.044 


13.057 


5 


13.380 


13.395 


13.300 


13.322 


6 


13.710 


13.466 


13.345 


13.322 


7 


13.810 


13.483 


13.552 


13.590 


8 


13.959 


13.717 


13.586 


13.590 


9 


14.448 


13.825 


13.680 


13.751 


10 


14.719 


14.042 


13.745 


13.751 


11 


- 


14.905 


14.985 


14.944 


12 


- 


15.105 


15.012 


15.058 


13 


- 


15.129 


15.057 


15.058 


14 


- 


15.147 


15.109 


15.172 


15 




15.287 


15.125 


15.172 


16 




15.381 


15.187 


15.172 


17 




15.439 


15.309 


15.322 


18 




15.451 


15.396 


15.322 


19 




15.496 


15.421 


15.435 


20 




15.496 


15.433 


15.435 



lattice spacing, f3 — T/h, the temperature is related to j3 via r = l/(Pks), 
and ks is the Boltzmann constant. 

Since we have approximated H by H e g, we can express those thermodynamical 
observables via the eigenvalues of the effective Hamiltonian 



N 
n=l 

— cfr N E eS e-P E ™ s 

E w= Jik~z s W' 
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Fig. 1. Average energy U of the Klein-Gordon model on a 1+1 dimensional lattice. 
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Fig. 2. Specific heat (C/ks) of the Klein-Gordon model on a 1+1 dimensional lattice. 



C cS (f3) = k B f3 2 



N (Effe-^ 



\n=l 



(14) 



Here N is the maximum number of eigenvalues available for a given N stoc h- 
Since this is a static system, the eigenvalues should in principle not vary with f3. 
This is observed numerically within statistical errors when f3 and N stoc h are not 
too small. Using this assumption and the spectrum at j3 — 2, we can compute 
thermodynamical quantities for other values of (3. For the N stoch = 1000 case, 
the average energy as a function of f3 is shown in Fig.[l] and the specific 
heat is shown in Fig. [2]. The results from the MC Hamiltonian are in good 
agreement with the analytical results for /3 > 1. However, for (3 < f3i ow ~ 1, 
which corresponds to the temperature r > r up = (kBPiow)" 1 , then the specific 
heat and also the average energy obtained from the MC Hamiltonian start 
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to deviate from the analytical result. This shows that the MC Hamiltonian is 
valid only in some finite temperature window. 

6 Summary 

In this paper, we have extended the Monte Carlo Hamiltonian method with a 
stochastic basis to quantum field theory (QFT), and taken the Klein-Gordon 
model as an example. The results are very encouraging. We believe that the 
application of the algorithm to more complicated systems will be very inter- 
esting. It will allow a non-perturbative investigation of physics beyond the 
ground state. 
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